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ABSTRACT 
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We present a model of y-ray emission through neutral pion production and decay in two- 
temperature accretion flows around supermassive black holes. We refine previous studies of 
such a hadronic y-ray emission by taking into account (1) relativistic effects in the photon 
transfer and (2) absorption of y-ray photons in the radiation field of the flow. We use a fully 
general relativistic description of both the radiative and hydrodynamic processes, which al- 
O . lows us to study the dependence on the black hole spin. The spin value strongly affects the 

y-ray emissivity within ~ 10 gravitational radii. The central regions of flows with the total 
^ I luminosities L< 10""' of the Eddington luminosity (LEdd) are mostly transparent to photons 

with energies below 10 GeV, permitting investigation of the effects of space-time metric. For 
such L, an observational upper limit on the y-ray (0.1 - 10 GeV) to X-ray (2-10 keV) lumi- 
nosity ratio of io.i-iOGev/i'2-iokeV ^0.1 can rule out rapid rotation of the black hole; on the 
^ ' other hand, a measurement of L().i-ioGev/i2-iokeV ~ 0.1 cannot be regarded as the evidence 

of rapid rotation, as such a ratio can also result from a flat radial profile of y-ray emissivity 

(which would occur for nonthermal acceleration of protons in the whole body of the flow). 

1^ ' At L> lO^^Lfidd, the y-ray emission from the innermost region is strongly absorbed and the 

• . observed y-rays do not carry information on the value of a. We note that if the X-ray emission 

^i ' observed in Centaurus A comes from an accretion flow, the hadronic y-ray emission from the 

2X. ', flow should contribute significantly to the MeV/Ge V emission observed from the core of this 

object, unless it contains a slowly rotating black hole and protons in the flow are thermal. 
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1 INTRODUCTION may be a signature allowing to test their fundamental nature. The 

model of M97 relied on a non-relativistic ADAF model and their 
Early investigations of black hole accretion flows indicated that ^ ^. . j u /-m o » i . /-^rvm i r. 

^ computations were improved by Oka & Manmoto (2003; hereafter 

tenuous flows can develop a two-temperature structure, with pro- 
ton temperature sufficient to produce a significant y-ray luminosity 
above 10 MeV through n" production (e.g. Dahlbacka, Chapline & 
Weaver 1974). The two-temperature structure is an essential feature 



OM03) who used a fully general relativistic (GR) model of the flow. 
However, both M97 and OM03 neglected the Doppler and gravita- 
tional shifts of energy as well as gravitational focusing and cap- 
turing by the black hole, which is a major deficiency because the 
of the optically-thin, advection dominated accretion flow (ADAF) . . . , , , ,,,,,,,,,,■ 

•^ -^ y-ray emission is produced very close to the black hole s horizon, 

model, which has been extensively studied and successfully applied ^^ .i i .i i i ^ j ^i ■ ^ i i ^- r 

^ J i-t- Furthermore, both works neglected the internal absorption ot y-ray 

to a variety of black hole systems (see, e.g., reviews in Yuan 2007, 



photons to pair creation, which effect should be important in more 
luminous systems. 



Narayan & McClintock 2008, Yuan & Narayan 2013) over the past 
two decades, following the work of Narayan & Yi (1994). Ma- 
hadevan, Narayan & Krolik (1997; hereafter M97) pointed out that ADAFs are supposed to power low-luminosity AGNs, like 

y-ray emission resulting from proton-proton collisions in ADAFs Fanaroff-Riley type I (FR I) radio galaxies or low-luminosity 

Seyfert galaxies, and a measurement, or even upper limits on their 

y-ray emission, may put interesting constraints on the properties 

* E-mail: niedzwiecki@uni.lodz.pl (AN), fgxie@shao.ac.cn (FGX), of the source of high-energy radiation in such objects. M97 and 

agastepnik82@gmail.com (AS) OM03 considered only the CGiJO/EGRET source in the direction 
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of the Galactic Center for such an analysis. Significant progress in 
exploration of the y-ray activity of AGNs which has been made 
after their works, thanks to the Fermi mission, motivates us to de- 
velop a more accurate model of the hadronic y-ray emission from 
ADAFs. Detections of y-ray emission from objects with misaligned 
jets (e.g. Abdo et al. 2010b) are most relevant for our study. Their 
y-ray radiation is usually explained as a jet emission; we show that 
emission from an accretion flow may be a reasonable alternative, 
at least in some FR Is. We focus on modelling of radiation in 100 
MeV - 10 GeV energy range, relevant for the Fermi-hAT measure- 
ments of the FR I radio galaxies (Abdo et al. 2010b) and over which 
the upper limits in Seyfert galaxies are derived (Ackermann et al. 
2012). 

The dependence of the y-ray luminosity on the black hole 
spin parameter makes a particularly interesting context for such an 
investigation. Already a rough estimate by Shapiro, Lightman & 
Eardley (1976) indicated a strong dependence of the y-ray lumi- 
nosity from a two-temperature flow on the spin of the black hole 
and, then, they suggested that this effect may serve as a means to 
measure the spin value (see also Eilek & Kafatos 1983 and Colpi, 
Maraschi & Treves 1986). OM03, who made GR calculations for 
the modern ADAF model, found a dramatic dependence of the y- 
ray luminosity on the spin value in models with thermal distribution 
of proton energies, however, they concluded that the dependence is 
weak if protons have a nonthermal distribution. In this work we ex- 
tend the analysis of this issue and clarify some related properties. 

We find global solutions of the hydrodynamical ADAF model, 
which follows Manmoto (2000), and use them to compute the y- 
ray emission. Similarly to M97 and OM03 we take into account 
emission resulting from thermal and nonthermal distribution of pro- 
ton energies; we use similar phenomenological models, with some 
modifications which allow to illustrate separately effects due to lo- 
cal distribution of proton energies and to radial profile of y-ray 
emissivity. We also use our recently developed model of global 
Comptonization (Niedzwiecki, Xie & Zdziarski 2012; hereafter 
N12, see also Xie et al. 2010) to compute the X-ray emission, which 
allows to investigate the internal absorption of y-ray photons to pair 
creation in the flow. 

In our computations we assume a rather weak magnetic field, 
with the magnetic pressure of I/lOth of the total pressure, sup- 
ported by results of the magnetohydrodynamic (MHD) simulations 
in which amplification of magnetic fields by the magneto-rotational 
instability typically saturates at such a ratio of the magnetic to the 
total pressure (e.g. Machida, Nakamura & Matsumoto 2004, Hirose 
et al. 2004, Hawley & Krolik 2001). We investigate the dependence 
on the poorly understood parameter in ADAF theoiy, 6, describing 
the fraction of the turbulent dissipation that directly heats electrons 
in the flow. We take into account only one value of the accretion 
rate, but the considered ranges of the spin and 6 parameters yield a 
rather large range of bolometric luminosities of ~ 10"* to 10"^ of 
the Eddington luminosity. In our paper we present both the spectra 
affected by yy absorption and those neglecting the absorption ef- 
fect; the latter may be easily scaled to smaller accretion rates, for 
which the yy absorption becomes unimportant. 



2 HOT FLOW MODEL 

We consider a black hole, characterised by its mass, M, and an- 
gular momentum, J, surrounded by a geometrically thick accre- 
tion flow with an accretion rate, M. We define the following di- 
mensionless parameters: r = R/Rg, a = J/(cRgM), m = MIMeaa, 



where M-^a = L^mIc^, Rg = GM/c^ is the gravitational radius and 
LEdd = 4nGMmj,c/cr-T is the Eddington luminosity. Most results 
presented in this work correspond to M = 2 x 10* Mg, in Fig[5^ we 
present also results for M = 2 x 10* Mq. We consider lii = 0.1 
and three values of the spin parameter, a = 0, 0.95 and 0.998. 
The inclination angle of the line of sight to the symmetry axis is 
given by Oobs- We assume that the density distribution is given by 
p(R,z) = p(/?, 0)exp(-z^/2//^), where H is the scale height at r. 
We assume the viscosity parameter of o- = 0.3 and the ratio of the 
gas pressure (electron and ion) to the total pressure of j8b = 0.9. 
The fraction of the dissipated energy which heats directly electrons 
is denoted by 6. 

Our calculations of hadronic processes are based on global so- 
lutions of the fully GR hydrodynamical model of two-temperature 
ADAFs, described in N12, which follows closely the model of 
Manmoto (2000). Here we recall only the ion energy equation, 
which is most important for the present study: 



o = (i-<5)ev,s + eco, 



A,e - Sir 



(1) 



where Aje is the Coulomb rate, the compressive heating and the 
advection of the internal energy of ions, respectively, are given by 



Geo, 



Mpi dlnp 



e,n 



Mp, 



dlnTi 



'™"'" iTTRp dR ' "^"" 27TRp{r,-l) dR 
and the viscous dissipation rate, per unit area, is given by 



e,i, = -apH(2n) 



r' dr ' 



(2) 



(3) 



where p = (pi + Pc)IP^, P\ is the ion pressure, p^ is the electron 
pressure, Fj is the ion adiabatic index, Q is the angular velocity 
of the flow, y^ is the Lorentz factor of the azimuthal motion and 
A = r^ + r^c? + 2ra. The form of the energy equation given in equa- 
tion ([T) is standard in ADAFs theory, although actually it should 
include an additional term describing the direct cooling of protons 
to pion production, Qy. In our calculation of hadronic processes we 
find that Qy is approximately equal to Aje at r < 10. At m = 0.1, 
considered in this work, both Qy and Aic are much smaller, by over 
3 orders of magnitude, than the effective heating Q^\^ + 2compr and 
the heating is fully balanced by the advective term, 2i„,. This justi- 
fies our neglect of the direct hadronic cooling. 

The only difference between our GR model and that 
of Manmoto (2000) involves the simplifying assumption of 
d \n{R)ld ln{H) = 1 adopted in the latter; we do not follow this 
simplification and an exact H{R) profile is considered in all our 
hydrodynamic equations. We note that the simplification has a con- 
siderable effect in the central part of the flow, e.g. it results in an un- 
derestimation of the proton temperature by a factor of ~ 1.5 within 
the innermost \QR^. Applying the above simplifying assumption 
we get exactly the same flow parameters as Manmoto (2000); note, 
however, that Manmoto (2000) assumed an equipartition between 
the gas and plasma pressures, withjSB = 0.5, which in general gives 
a smaller proton temperature than /Jb = 0.9 assumed here. In par- 
ticular, for a = and 6 = 10"', models with /3b = 0.9 give the 
proton temperature larger by a factor of ^ 4, close to the horizon, 
than models with/?B = 0.5. This underlies also the differences in 
the y-ray luminosity levels between the thermal models of Oka & 
Manmoto (2003) and ours, as discussed in Section 3. 

To obtain global transonic solutions we have to adjust the spe- 
cific angular momentum per unit mass accreted by the black hole, 
for which the accretion flow passes smoothly through the sonic 
point, Ts. We note that this condition permits for two kinds of so- 
lutions, below referred to as a 'standard' and a 'superhot' solution. 
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Figure 1. Radial profiles of the dissipative heating rates, Qvis (a), the proton temperature (b) and the proton number density (c) of our hot-flow solutions for 
(5 = 10" . In all panels, the dashed (red) lines are for a = 0.998, the solid (black) lines are for a = 0.95 and the dotted (blue) lines are for a = 0. In panel (a), 
Qvis denotes a vertically integrated rate, so 2vis^^ gives the heating rate (per unit volume) times volume. The green (dashed) line in panel (a) shows Qyjs in 
the superhot solution (see text) for a = 0.998. M = 2 X 10 Mq, m = 0.1, a = 0.3 and /3b = 0.9 in this and all further figures in this paper. 



The latter (superhot) has much larger proton temperature and den- 
sity, furthermore, the sound speed is large and the sonic point lo- 
cated in the immediate vicinity of the event horizon, e.g. r, =: 1.2 
for a = 0.998. In the standard solutions the sonic point is located at 
larger distances, r, > 2. Taking into account rather extreme proper- 
ties of the superhot solutions (specifically, a very large magnitude 
of 2vis illustrated in Fig.[TJa) and discussed in Section lOl we ne- 
glect them in this work and for all values of a we consider only the 
standard solutions which are consistent with solutions of the model 
investigated in several previous studies (e.g. Manmoto 2000, Yuan 
et al 2009, Li et al. 2009). Note, however, that in our previous works 
(N12, Niedzwiecki, Xie & Beckmann 2012) we considered the su- 
perhot solution with a = 0.998, then, the results for a = 0.998 
discussed in those works correspond to flows with larger proton 
temperature and density (both by a factor of ~ 5) than these con- 
sidered in the present study. 

Fig. [T] shows the dependence on the black hole spin of some 
parameters from our solutions which are crucial for the hadronic 
y-ray production. Rotation of the black hole stabilizes the circular 
motion of the flow which yields a higher density (through the conti- 
nuity equation). Furthermore, the stabilized rotation of the flow re- 
sults in a stronger dissipative heating giving a larger proton temper- 
ature for larger a. All these differences are significant only within 
the innermost ~ lORr,. 



3 HADRONIC y-RAY EMISSION AND RELATIVISTIC 
TRANSFER EFFECTS 

The hydrodynamical solutions set the proton density, iij,, and tem- 
perature, Tp, as a function of radius. In principle, this should allow 
to determine the y-ray emissivity, resulting from neutral pion pro- 
duction in proton-proton collisions and their subsequent decay into 
y-ray photons, in the rest frame of the flow. However, details of 
this process are subject to an uncertainty related to the distribution 
of proton energies, which is unlikely to be thermal in optically thin 
flows (see discussion in Section |6^ . Following M97 and OM03 we 
assume that the temperature from the global solution functions as 
a parameter specifying the average energy of protons in the plasma 
which, however, does not have to have a thermal distribution. We 
consider several phenomenological models which must satisfy the 
obvious requirements that at each radius (1) the number density of 



protons equals «p(r), determined by the global ADAF solution and 
(2) the average energy of protons equals the average energy 



f/thC^p) = e„m„c\6 + 150p)/(4 + 5e„). 



(4) 



of the Maxwellian proton gas with temperature, Tj,(r), determined 
by the global ADAF solution, where dp = kTp/mpC' and we use 
the simplified (cf. Gammie & Popham 1998) relativistic form of 

f/.h(ep). 

We consider models involving various combinations of ther- 
mal 



n,h(y) = «,hy'yeexp(-y/(?p)/[(?p/f2(l/ep)], 
and power-law 
nni(y) = n„,is- l)y"', 



(5) 



(6) 



distributions of proton energies, where n,h and Wpi are the local den- 
sities of these two populations. The thermal model (model T) as- 
suming a purely Maxwellian distribution of protons and the non- 
thermal model (model N, same to nonthermal models of M97 and 
OM03), assuming that the total energy is stored in the power-law 
distribution of a small fraction of protons, allow us to estimate the 
minimum and maximum level of y-ray luminosity, respectively, for 
a given set of {M, m, a, a, /}, 6). Mahadevan (1999) and OM03 con- 
sidered the model involving the mixture of the thermal and power- 
law distributions, with the radius-independent parameter character- 
izing the fraction of energy that goes into the two distributions. 
Deviations of such a model from model N are trivial, with the y- 
ray luminosity linearly proportional to the fraction of energy going 
to the power-law distribution. In this work we consider a different 
hybrid model (model H) with the radius dependent normalization 
between the power-law and the thermal distribution, which allows 
us to illustrate some additional effects. 

The detailed assumptions on the parameters of these models 
are as follows («p(r) and Tj,(r) denote values given by the global 
ADAF solution): 

Model T assumes a purely Maxwellian distribution of protons, 
equation l|5j, with n,h = Wp(r) and dp = kTp(r)/mpC^. 
Model N assumes that a fraction t// of protons form the power- 
law distribution, equation (|6}, with the radius-independent index 
.? and rip] = i/f(r]np(r), and the remaining protons are cold, with the 
Lorentz factor y ^ 1 (y = 1 is assumed in the computations). The 
radius-dependent fraction ^ is determined by 
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Figure 2. Dashed (green) lines show the radial profiles of the vertically-integrated y-ray emissivities, Qy, for models with a 
show the local contribution to the luminosity at infinity from a unit area of the flow neglecting the yy absorption, /iunabs> for i^ 
the local contribution to the luminosity at infinity from a unit area of the flow taking into account the yy absorption, /labs, for o = 
and (lower, blue), (a) model T; (b) model H with .v = 2.6; (c) model N with s = 2.6. All models assume S = 10"'. 



= 0.95. Sohd (magenta) lines 
= 0.95. The dotted lines show 
0.95 (upper, black) and = 



i/rmpC 



(/,h [T,(r)] . 



(7) 



Model H assumes that an efficient nonthermal acceleration oper- 
ates only within the central ~ 15R^, where the average proton ener- 
gies resulting from the ADAF solutions become relativistic. Specif- 
ically, we assume that at each radius at r < 15 a fraction of pro- 
tons form a thermal distribution at a subrelativistic temperature of 
T = 4.3 X 10" K (dp = 0.04), and the remaining form a power-law 
distribution (equation |6ll with a constant (i.e. radius-independent) 
index s and npi = tf/(r)np(r). The relative normalization of these 
two distributions is determined by 



il/m„c^ 6.6 T r n 

Try + ^1 - ^)^'"p^ = f^.h [^p('-)] 



(8) 



(where the factor 6.6/4.2 results from equation (|4) with dp = 0.04). 
At r > 15, where T^ < 4.3 x 10" K, there are no non-thermal 
protons in this model, which then results in a negligible pion pro- 
duction at such distances, similar as in model T. The chosen value 
of r = 4.3 X 10" K gives a smooth transition between a purely ther- 
mal and a hybrid plasma at r = 15, however, radiative properties of 
model H are roughly independent of the specific value of the tem- 
perature of the subrelativistic thermal component. We remark also 
that T = 4.3 X 10" K is close to the limiting temperature above 
which the pion production prevents thermalization of protons (see 
Stepney 1983, Dermer 1986b) 

The efficiency of pion production by protons with the power- 
law distribution increases with the decrease of the power-law index 
.?. On the other hand, the fraction if/ decreases with decreasing s, 
roughly asifi oc (.s-2). These two effects balance each other yielding 
the largest luminosity in 0. 1-10 GeV range, Lo.i_iogcV> f^or s ^ 2.5- 
2.6. For 2.3 < .? < 2.8, the dependence of Lq.i-iogcV on i- is weak; 
for i- = 2.1, Lo.i-ioGeV is by a factor of ~ 2 smaller than for s = 2.6. 
To estimate the maximum value of Lo.i_iogcV that can be produced 
in a flow with given parameters, in our computations for models N 
and H we use s = 2.6. For all values of a, dp > 0.1, and also i/r > 0. 1 
in models H and N with s = 2.6, within the innermost several R„. 

In our solutions of the flow structure we assume that protons 
are thermal and we use the thermal form of the gas pressure. Then, 
our models N and H with non-thermal proton distributions are not 
strictly self-consistent, as their pressure may deviate from the ther- 
mal prescription. However, this is a rather small effect, e.g. the 
pressure of the purely non-thermal distribution (model N) differs 
by 20-30 per cent from the pressure of a thermal gas with the same 
internal energy. 



For a given distribution of proton energies we determine the 
y-ray spectra in the flow rest frame, strictly following Dermer 
(1986a, 1986b), in a manner similar to M97 and OM03; however, 
we do not apply the following simplification underlying their non- 
thermal model. As argued in M97, the fraction of nonthermal pro- 
tons should be small, i/' <c 1, and, therefore, interactions of non- 
thermal protons with other nonthermal protons may be neglected; 
hence, only interaction of nonthermal protons with cold protons are 
taken into account in their computations. We remark that such an 
approach underestimates the y-ray luminosity, e.g. by a factor of 
~ 2 in model N with a = 0.95 and s = 2.6 (for which if/ ^ 0.4 
in the innermost region). In all our models we take into account 
interaction of protons with all other protons. 

To compute the y-ray luminosity and spectra received by dis- 
tant observers we use a Monte Carlo method similar to that de- 
scribed in Niedzwiecki (2005). We generate y-ray photons isotrop- 
ically in the plasma frame, make a Lorentz transformation from the 
flow rest frame to the locally non rotating (LNR) frame and then 
we compute the transfer of y-ray photons in curved space-time; 
see, e.g., Bardeen et al. (1972) for the definition of LNR frames 
and the equations of motion in the Kerr metric. 

The dashed lines in Fig. [2] show the radial profiles of the 
vertically-integrated y-ray emissivity, Qy (Qy gives the energy 
emitted from the unit area per unit time) for models T, H and N 
with a = 0.95. The solid lines in Fig. |2] show the radial profiles of 
the vertically-integrated local luminosity (the energy per unit time 
reaching infinity from the unit area at a given r). The local lumi- 
nosity profiles shown by the solid lines neglect the yy absorption, 
so the difference between the dashed and solid lines is only due 
to the relativistic transfer effects. Fig. [3] shows the coiTesponding 
y-ray spectra and compares them with the spectra for a = 0. At 
r < 10 both models N and H are characterised by similar values of 
if/ and produce similar amounts of y-rays. In both models T and H 
the contribution from r > 10 is very weak; in model N the radial 
emissivity is much flatter despite if/ being small, e.g. if/ < 5 x 10"' 
at r > 100. Comparing models T and H we can see the effect of the 
local proton distribution function and by comparing models H and 
N we can see the effect of the radial emissivity. 

For the thermal distribution of protons, the rest-frame photon 
spectra are symmetrical, in the logarithmic scale, around ~ 70 MeV 
but in EFe units they peak around 200 MeV; the position of the 
maximum in the spectra observed by distant observes is slightly 
redshifted. Note that the difference of y-ray luminosities, Ly, be- 
tween a = and 0.95 in our model T is much smaller than that 
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Figure 3. Dashed (blue) and solid (black) lines show the rest frame and the 
observed y-ray spectra, respectively, for model T (ab) and model H with 
s = 2.6 (cd). Dotted (red) and dot-dashed (magenta) lines in (cd) show the 
rest frame and the observed y-ray spectra, respectively, for model N with 
s = 2.6. All spectra are for 6 = 10"''; panels (a) and (c) are for a = 0, panels 
(b) and (d) are for a = 0.95. In this figure, the observed spectra neglect yy 
absorption, so they are effected only by GR effects. 
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Figure 4. Observation-angle dependent y-ray spectra taking into account 
and neglecting the yy absorption for model H with s = 2.6; a = 0.998 
(a) and a = (b). The dot-dashed (red) and dashed (blue) lines show the 
spectra observed at ©oj,, = 70°, with and without absorption, and the solid 
(magenta) and dotted (black) lines show the spectra observed at 6^^,^ = 30°, 
with and without absorption, respectively. 



this work are quasi-spherical and optically thin and hence their ap- 
pearance depends on the viewing angle primarily due to the rel- 
ativistic transfer effects. Most importantly, trajectories of photons 
emitted close to a rapidly rotating black hole are bent toward its 
equatorial plane. Therefore, the y-ray radiation has a significant in- 
trinsic anisotropy in models with large a, with edge-on directions 
corresponding to larger y-ray fluxes. 



derived by OM03, whose thermal models with a = and 0.95 give 
Ly differing by approximately three orders of magnitude. The dif- 
ference is due to different values of /Jg assumed here and by OM03, 
which result in different 6p, as discussed in Section |2] The depen- 
dence of Ly on dp changes around 9^ x 0.1 (see, e.g., fig. 3 in 
Dermer 1986b). At lower temperatures, Ly is extremely sensitive to 
Op, with the increase of Op by a factor of 2 yielding the increase of 
Ly by over two orders of magnitude. At 0p > 0.1, the dependence 
is more modest, e.g. the increase of Op from 0.2 to 0.4 results in 
the increase of Ly by only a factor of ~ 2. For ySa = 0.9 assumed 
in this work, 0p > 0.1 at small r for all values of a, making the 
y-ray luminosity much less dependent on the black hole spin. For 
j8b = 0.5, assumed by OM03, the proton temperature is small, with 
the maximum value of Op « 0.03 for a = 0, which leads to the 
strong dependence of Ly on a. 

For both model H and N, the spectrum at ii > 1 GeV has the 
same slope as the power-law distribution of proton energies. For 
model H with s = 2.6, Ly is by a factor of 3 larger than in model T. 
Rather small difference between Ly in our thermal and nonthermal 
models is again due to our assumption of a weak magnetic field. 
At smaller /3b, resulting in smaller Op, the presence of even a small 
fraction of non-thermal electrons leads to the increase of Ly by or- 
ders of magnitude, as can be seen by comparing the emissivities of 
our models N and T at r > 10 (see also M97). 

For models T and H the bulk of the y-ray emission comes from 
r < 10 (Fig. [2^b) and the GR transfer effects reduce the detected 
y-ray flux by approximately an order of magnitude. In model N 
the magnitude of the GR effects on the total flux is reduced due 
to strong contribution from r > 10 (which is weakly affected by 
GR). Also in model N, the contribution from r > 10, which for 
a = approximately equals the contribution from r < 10, reduces 
the difference between the y-ray fluxes observed for a = and 
a = 0.95 to only a factor of ~ 2, see Fig.[3lcd). 

The viewing-angle dependent spectra for model H, which 
would be observed (if unabsorbed) by distant observers, are shown 
by the dashed and dotted lines in Fig. |4] The flows considered in 



4 COMPTONIZATION AND yy ABSORPTION 

The absorption of y-rays in the radiation field of the flow has been 
calculated in a fully GR model by Li et al. (2009). Here we use a 
similar approach with the major difference involving the computa- 
tion of target photon density. Li et al. (2009) considered the prop- 
agation of photons with energies of 10 TeV, which are absorbed 
mostly in interactions with infra-red photons. Those low energy 
photons are produced primarily by the synchrotron emission which 
can be simply modelled using its local emissivity. In turn, pho- 
tons with energies in 0.1-10 GeV range, considered in this work, 
are mostly absorbed by the UV and soft X-ray photons, which are 
produced by Comptonization. Then, an exact computation of the 
angular-, energy- and location-dependent distribution of the target 
photon field requires the precise modelling of the Comptonization 
taking into account its global nature. In our model we apply the 
Monte Carlo (MC) method, described in detail in N12, with seed 
photons for Comptonization from synchrotron and bremsstrahlung 
emission. 

We find self-consistent electron temperature distributions us- 
ing the procedure described in N12; we iterate between the solu- 
tions of the electron energy equation (analogous to equations [T}{3j 
note that here we include the direct electron heating, while N12 as- 
sumes 6 = 0) and the GR MC Coinptonization simulations until 
we find mutually consistent solutions. In Fig.[5]we show the result- 
ing spectra. Fig. |6] shows the radial profiles of the radiative cooling 
of electrons (strongly dominated by Comptonization), i2compt> for 
6 = 10"'' and compares them with the y-ray emissivity, Qy, for 
model T. Note that Qy is much steeper than Qcompt so the GR ef- 
fects are more important for the y-ray than for the X-ray emission. 

As we can see in Fig. |6] and also in the corresponding spec- 
tra in Fig. |5^, for 6 = lO"'' the black hole spin negligibly affects 
the Comptonized radiation; this property results from a large mag- 
nitude of the compression work, which is roughly independent of 
a and dominates the heating of electrons for small values of d (cf. 
N12). 
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Figure 6. Radial profiles of the y-ray emissivity, Qy (for model T), and 
the Comptonization rate Qcompt> for « = 0.998 (dashed, red), 0.95 (solid. 
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Figure 5. Angle-averaged spectra received by a distant observer; the synchrotron and Comptonized (radio to X-rays) and hadronic (y-rays) components are 
shown separately. In both panels a = 0.998 (dashed, red), 0.95 (solid, black) and (dotted, blue) (a) Models with S = 10"^; the y-ray spectra con'espond to 
model H with s = 2.6. (b) Models with S = 0.5; the y-ray spectra correspond to model N with s = 2.6 for a = and 0.998 and to model H with s = 2.6 for 
a = 0.95. The lower pair of lines in panel (a), rescaled by a factor of 10, show the spectra (leptonic component) for M = 2 X lO^M©. 



For (5 > 0.1 the direct heating contributes significantly to the 
heating of electrons and for S = 0.5 it dominates over other heating 
processes at r < 100 for all values of a. Then, the dependence of 
2vis on a results in a noticeable dependence of the Comptonized 
radiation on a for 5 > 0. 1 (see also Xie & Yuan 2012 for a recent 
study of the dependence of X-ray luminosity on 6). The radiative 
efficiency increases from 77 = 0.004 for all values of a at 5 = 10"' 
tori = 0.02 foia = 0,r] = 0.08 for a = 0.95 and r] = 0.1 for a = 
0.998 at 6 = 0.5. Despite considering only one value of accretion 
rate, our solutions span a range of bolometric luminosities, from 
L a 4 X lO^^'LEdd (for 6 = 10"') to L r; lO^-Lndd (for 6 = 0.5 and 
a = 0.998). The corresponding X-ray spectral slopes harden from 
Fx - 1.7 to Fx - 1.5 with the increase of L. Note that these values 
correspond to the range of parameters close to the turning point 
in the L-T correlation observed in AGNs (e.g. Gu & Cao 2009). 
Then, we likely consider here the range of the largest luminosities 
of the flows in which synchrotron emission is the dominant source 
of seed photons for Comptonization (see discussion and references 
inN12). 

Having found the self-consistent solutions, described above, 
we apply our MC model to tabulate the distribution of all pho- 
tons propagating in the central region (up to r^^t = 1000), 
dnj,t,(R,6,E]j^,Clijq)/dE]^f^dQ.]j^ (in photons cm"' eV"' sr"'), 
where R and 6 are the Boyer-Lindquist coordinates, £ln is the pho- 
ton energy in the LNR frame and d^LN is the solid angle element 
in the LNR frame. 

To compute the optical depth to pair creation, Tyy, we closely 
follow the method for determining an optical depth to Comp- 
ton scattering in the Kerr metric, see Niedzwiecki (2005) and 
Niedzwiecki & Zdziarski (2006), however, here we calculate the 
probability of pair creation in the LNR frame whereas for the 
Compton effect an additional boost to the flow rest frame is ap- 
plied. While Compton scattering is most conveniently described in 
the plasma rest frame, pair production can be simply modelled in 
the LNR frame and, thus, the transformation to the flow rest frame 
is not necessary here. We solve equations of the photon motion in 
the Kerr metric and we determine the increase of the optical depth 
along the photon trajectory from 



black) and (dotted, blue) in models with 6 - 
integrated rates. 



in the LNR frame and 6(ln is the angle between the interacting 
photons in the LNR frame. 

GR affects the 77 opacity through (1) bending the trajectories 
of both the 7-ray photon and target photons and (2) changing en- 
ergies of both the y-ray photon and target photons. As an example, 
the neglect of the gravitational shift of the y-ray photon energy, by 
using (ryy(E]^f^,Ey, (Pln) (where Ey is the energy at infinity) instead 
of cryy(£LN> Eyi^fi, (^LN ) 111 equatlou l|9), underestimates r-^y by a fac- 
tor of X 2-3 for photons emitted from the innermost region. 

In Fig. [7] we show values of the total optical depth, Tyy(r), in- 
tegrated along the outward radial direction in the equatorial plane 
from the emission point at the radial coordinate r to the outer 
boundary at ro^f As we can see, the 77 opacity is a strong func- 
tion of both the y-ray energy and the location in the flow. The dot- 
ted lines in Fig. [2] show how the 77 absorption attenuates 7-rays 
observed from a given r. 

It is apparent that around ra ~ 0. 1 flows undergo transition 
from being fully transparent to mostly opaque to y-rays. In our 
models with the Eddington ratio L/Lem = 4 x lO"'*, the flow is 
fully transparent to photons with energies < 100 MeV; at higher 
energies the absorption leads to moderate attenuation, with the in- 
crease of the photon index at £ > 1 GeV by AF == 0.2, see Fig.|4] 
The size of the y-ray photosphere (the surface of Tyy = 1) increases 
with increasing L and for L =^ 10"^LEdd the GeV photons cannot 
Gould & Schreder 1967), £yLN is the energy of the y-ray photon escape from r < 10. At such L, our model H gives spectra with 
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where d/LN is the length element in the LNR frame, 
(Tyy(Ei^t^,EyLN,S]^N) Is thc pair production cross section (e.g.. 
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Figure 8. Observed, angle-averaged y-ray luminosity in 0.1-10 GeV range in model T (a) and model N with .s = 2.6 (b) and X-ray luminosity in the 2-10 keV 
range (c) as a function of 5 for a = 0.998 (circles), a = 0.95 (tiiangles) and a = (squares). The filled symbols in (ab) show Lo.l-lOGeV after the yy absorption, 
the open symbols show the luminosity Lunabs neglecting the absorption. 
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Figure 7. The optical depth to pair creation for radially outgoing y-ray pho- 
tons with Ey = 100 MeV, 1 GeV and 10 GeV from bottom to top, as a 
function of the radial distance of their point of emission for a = 0.95 and 
6 = 10"' are shown by the solid (black) lines. The dotted (blue) and dashed 
(red) lines are for a = 0.95 and in models 6 = 0.5; in these models Tyy is 
shown only for £" = 10 GeV for claiity. 



a clear cut-ofF around 1 GeV (see the solid line in Fig. [Sj)) which 
could be measured by Fermi. In other cases absorption leads to a 
smooth softening of the spectra. 

In terms of the 2-10 keV luminosity, La-iokcVi flows with 
i2-iokcV < lO^^LEdd should be fully transparent to MeV/GeV pho- 
tons. Flows with L2-iokeV > lO^^Lndd can emit significant amounts 
of unabsorbed y-rays only if their y-ray emissivities are strong at 
large r. E.g. in our model N, the luminosity of the flow at r > 50, 
which region would be outside the photosphere of 1 GeV photons 
even at much larger L2-iokcV ~ 10""LEdd> is ^oi-ioGcV - lO*" erg/s. 
Then, the y-ray luminosity exceeding lO*' erg/s can be expected at 
m > 0.3 if the y-ray emitting flow extends out to several tens of Rg, 
which propeity is, however, unclear as objects with high luminosi- 
ties often show signs of a cold disc extending to rather small radii 
(so the transition between the hot and cold flow may occur within 
the y-ray photosphere). Note that for such a scenario, with y-ray 
emission from a hot flow at large L, we expect a small luminosity 
ratio of Lo.i-ioGev/i2-iokev ~ 10"' regardless of the value of a. 



5 X-RAY VS y-RAY LUMINOSITY 

In Fig. [8] we summarize our results regarding the relation between 
the X-ray and y-ray luminosities. The range of expected Lo.i_iogcV 
is constrained from below by values indicated in Fig.|8}i for model 
T, and from above by values in Fig.[8j3 for model N. As we can see, 
the models give the luminosity ratios io.i-ioGev/i^-iokeV between 
~ 0.002 and 0.2. 

In model T, Lq i_iogcV for a = and a = 0.998 differ by a fac- 
tor of several; the unabsorbed luminosities difi"er by over an order of 
magnitude but for L close to 10"^LEdd the yy absorption reduces the 
difference to a factor of ~ 4. In model N, Lq.i-iogcV for a = and 
a = 0.998 differ by only a factor of ~ 2; as noted before, the differ- 
ence is reduced here due to contribution from large r. Model N for 
a = gives similar Lq.i-iogcV as model T with large a; larger den- 
sity and average energy for large a is approximately compensated 
by a larger fraction of protons above the pion production threshold 
for model N. Note, however, that - despite similar luminosities - the 
spectra for these two regimes differ significantly, see Fig.|3] 

In model H, Lo i-iogcV has a similar dependence on a as in 
model T, with a large difference between small and high values of 
a. We conclude that it is the radial distribution of y-ray emissivity, 
rather than the local proton distribution function, which reduces the 
dependence on a. This is even more hindering for attempts of as- 
sessing the spin value basing on the y-ray luminosity level. The pro- 
ton energy distribution function is reflected in the produced spectral 
shape and, therefore, it could be constrained observationally and 
taken into account in this kind of analysis. On the other hand, one 
cannot expect to derive information on the radial emissivity pro- 
file from observations, so an investigation of black hole spin values 
using the y-ray luminosity should be subject to significant uncer- 
tainty. 

In models with a dominant contribution from the central ~ 
lORg, intrinsic y-ray luminosities of flows around submaximally 
(a = 0.95) and maximally (a = 0.998) rotating black holes differ 
by a factor of ~ 2. For the latter, the flow extends to smaller radii 
and hence a larger proton temperature and density are achieved, see 
Fig.[T] However, a strong contribution from r > 10 as well as yy 
absorption make the difference between the observed luminosities 
insignificant. 

Regarding the dependence on 6, we notice a somewhat sur- 
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prising property of models with 6 = 0.5, which predict a larger 
y-ray luminosity than models with smaller 6 (the effect is more 
pronounced for larger values of a). The physical reason is that for 
6 = 0.5, the slight decrease of proton temperature is outweighed by 
the increase of density (through the decrease of both the radial ve- 
locity and the scale height with decreasing temperature), cf. Man- 
moto (2000), which leads to the increase of Ly oc jJpTp. At smaller 
values of d, the intrinsic L-y depends negligibly on 6, in agreement 
with the results of OM03 for a = and (5 < 0.3 in their nonthermal 
model. 

The scaling of density with ih and M in our global GR solu- 
tions only weakly differs from that of self-similar ADAF model, 
i.e. n oc m/M (e.g. Mahadevan 1997). Then, the intrinsic y-ray lu- 
minosity can be estimated as Lo.i-ioGeV - ^uiiabs(™/0.1)^(M/2 x 
10** Mq)"' erg/s, where L^nai,, is the unabsorbed luminosity shown 
by the open symbols in Fig.[8jab). This gives also the observed lu- 
minosity at m < 0.1, when the absorption effects are unimportant. 
Obviously, the above scaling with m is not relevant for m>0. 1, 
for which the increase of lii results in the decrease of the observed 
io.i-ioGev due to yy absorption. On the other hand, the linear scal- 
ing of the y-ray luminosity with M holds even when the yy absorp- 
tion is important. The spectral distribution of Comptonized radia- 
tion changes slightly with the black hole mass (due to the change 
of the synchrotron emission, which affects the position of Comp- 
tonization bumps, see Fig 5a), however, the effect is insignificant 
for the yy opacity and we have checked that it negligibly affects 
the observed y-ray spectra. 



6 DISCUSSION 

The form of the dissipation rate and the proton distribution func- 
tion are two major uncertainties for predicting the y-luminosity. 
We briefly discuss here the related effects in accretion flows. 



6.1 Viscous heating 

The form of the viscous dissipation rate given by equation ([!) re- 
sults from the usual assumption that the viscous stress is propor- 
tional to the total pressure, with the proportionality coefficient a. 
We use this form of the viscous stress for computational simplic- 
ity, however, as we discuss in N12, this may be an oversimplified 
approach (although it has support in MHD simulations, see below). 
Below we briefly compare the dissipation rate in our models with 
that predicted by the classical Novikov & Thome (1973) model. 
We define 2,;, integrated over the whole body of the flow as the 
total dissipation rate, 2vis,toti in our solutions and compare it to the 
dissipation rate in a Keplerian disc for the corresponding value of a, 
2NT,tot> given by the Novikov & Thorne model. Obviously, 2vis,tot 
does not need to match 2NT,tot closely. The latter value is calculated 
under the assumption that the shear stress vanishes at the radius 
of the innermost stable circular orbit (ISCO), which condition is 
not applicable to geometrically thick ADAFs and its release should 
lead to stronger dissipation. On the other hand, ADAFs are sub- 
Keplerian which property decreases the dissipation rate. 

Our 2vis,iot does not differ significantly from 2NT,tot! however, 
the efficiency of dissipation in our model comparatively increases 
with increasing a. Specifically, Qvis.tot = 0.025Mc^ ~ 0.52NT.tot 
for a = 0, !2vis,tot = 0.2Mc^ a Qmm for a = 0.95 and gvis.tot = 
0.58Mc^ « 1.52NT,iot for a = 0.998. Note that in an ADAF with 
a large value of a most of the dissipation occurs very deep in the 



potential, where relativistic effects strongly reduce the energy es- 
caping to infinity, whereas in a Keplerian disc the dissipation is 
less centrally concentrated and its radiation is subject to less severe 
reduction. For example, the gravitational redshift alone (neglecting, 
e.g., the photon capture under the event horizon) would give similar 
luminosities received by distant observers (^ 0.3Mc^) both in our 
ADAF and in Novikov & Thome models with a = 0.998, if all of 
the dissipated energy were converted into radiation. Obviously, in 
optically thin ADAFs only a small part of the dissipated energy is 
radiated away and most of it is accreted by the black hole, so the 
radiative efficiency is much smaller than 0.3. 

As we mentioned in Section |2] the model formally allows for 
two solutions and the above values of Qvis.tot correspond to the stan- 
dard solution, considered in previous sections. The superhot solu- 
tion has an extreme dissipation, with (2vis,toi - 5-lOMc-. This is 
not necessarily an unphysical property (note that gvis.iot is defined 
in the rest frame of the flow), as the observed luminosity of such 
flows does not exceed the accretion rate of the rest mass energy. 
Nevertheless, the magnitude of Qvis suggests that the underlying 
assumption of Qvis '^ P breaks down at very large p. 

The issue of the proper description of the Qvis term could be 
partially resolved by comparing analytic models such as our, aim- 
ing at a precise calculation of the produced radiation, with MHD 
simulations. The latter currently neglect radiative cooling, or use 
very approximate descriptions for it, in turn, they provide more 
accurate accounts of the dissipation physics. Such MHD simula- 
tions support some properties of our model, e.g. the Qvis '^ P pre- 
scription of viscous heating (Ohsuga et al. 2009). Furthermore, the 
GR MHD simulations have shown that flows cross the ISCO with- 
out any evidence that the shear stress goes to zero, which leads 
to the increase of the radiative efficiency, and the deviations from 
Novikov & Thome model increase with increasing H/R ratio (e.g. 
Noble, Krolik & Hawley 2009, Penna et al. 2010). We could not, 
however, quantitatively compare our models with such simulations, 
as the published results have a much smaller aspect ratio than our 
solutions {H/R > 0.5). 

Lastly, we remark that Gammie & Popham (1998) and 
Popham & Gammie (1998) present a model similar to ours but with 
a more elaborate description of the shear stress (in turn, they neglect 
radiative processes). We note that their results indicate a similar in 
magnitude stabilizing effect of the rotation of black hole. 

6.2 Proton distribution function 

As pointed out by Mahadevan & Quataert (1998), Coulomb colli- 
sions are too inefficient to thermalize protons in optically thin flows 
and the proton distribution function is determined by the viscous 
heating mechanism, which is poorly understood. Our solutions, 
with M = 2 X 10^ Mq and m = 0. 1, give the accretion time-scale, 
?a, much shorter than the proton relaxation time-scale, ?pp; e.g. at 
r < 20, fa < 10 hours and ?pp > 10^ hours. Clearly, the protons 
cannot redistribute their energy through Coulomb collisions. 

In solar flares, the best observationally studied example of 
particle acceleration/heating in a magnetised plasma, a significant 
fraction of the released energy is carried by non-thermal, high en- 
ergy particles (e.g. Aschwanden 2002), which strongly motivates 
for considering the nonthermal distribution of protons, as origi- 
nally proposed in M97. Applying the generic description of particle 
acceleration, see e.g. section 3 in Zdziarski, Malzac & Bednarek 
(2009), we check whether the conditions in ADAFs allow for pro- 
ton acceleration to ultrarelativistic energies, as assumed in our 
computations for the power-law distributions. The magnetic field 
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strength in our ADAF solutions is fi ^ 10 G, 100 G and 1000 G 
at r = 100, 10 and 2, respectively. Assuming an acceleration rate 
dE/dt oc ^eB, where e is the elementary charge and ^ is the accel- 
eration efficiency (^ < 1), we find that the maximum Lorentz factor 
limited by the synchrotron energy loss is y^^,, ~ 10* - 10^, depend- 
ing on r. Another condition, of the Larmor radius being smaller 
than the acceleration site size, i^acc, gives 7„,ax ~ 10' at r » 100, 
and larger values of y^.^,^ at smaller r, even if we safely assume 
^acc = l^g (= 3 X 10'^ cm). We conclude that the central region of 
a hot accretion flow may be a site of the acceleration of protons to 
energies which easily allow hadronic emission of photons even in 
the TeV range. This conclusion remains valid for the whole relevant 
range of accretion rates and masses of supermassive black holes, as 
the strength of the magnetic field scales as fi oc m'^^M"''^. 

Considering processes which could compete with proton- 
proton interactions, we note that proton-photon interactions are 
much less efficient in the central region of the flow. Namely, the 
number density of y-ray photons, which can elfectively interact 
with protons in photomeson production, is by a factor of ~ 10^ 
smaller than the number density of protons within the innermost 
few R„; the number density of hard X-ray photons, which can in- 
teract with protons in photopair production, is similar to «p. The 
cross-section for both channels of proton-photon interactions is by 
over two orders of magnitude smaller than the cross-section for 
proton-proton interaction, making the photo-hadronic production 
of secondary particles negligible. 



7 COMPARISON WITH OBSERVATIONS 

We briefly compare here predictions of our model with y-ray ob- 
servations of objects which may be powered by ADAFs. We note, 
however, that for a more detailed comparison additional physi- 
cal processes, related to charged pion production (cf. Mahadevan 
1999), should be taken into account. In particular, relativistic elec- 
trons produced by the pion decay should be important in modelling 
the emission in the MeV range. We also tentatively discuss the ori- 
gin in the accretion flow of the very high energy radiation detected 
from M87 and Sgr A*, for the latter under assumption that the con- 
tributions from two separate sources dominate above 100 GeV (Sgr 
A*) and at lower energies (diffuse emission) in the radiation ob- 
served from the Galactic Center region. However, we note that the 
opacity at very high energies may be affected by the nonthermal 
synchrotron emission of the relativistic electrons. The work on the 
model implementing the effects of charged pions is currently in 
progress. 



7.1 Misaligned AGNs 

The misaligned AGNs detected by the Fermi-hKT include seven PR 
I radio galaxies and four FR lis (Abdo et al. 2010b). The low-power 
FR I galaxies are supposed to be powered by radiatively inefficient 
accretion flows (e.g. Balmaverde, Baldi & Capetti 2008) and, thus, 
are more relevant for application of our results. Interestingly, Wu, 
Cao & Wang (201 1) assess that supermassive black holes in FR Is 
rotate rapidly, with a > 0.9. The X-ray emission of more luminous 
FR Is is supposed to come from an accretion flow (e.g. Wu, Yuan 
& Cao 2007), but their y-ray emission tends to be interpreted in 
teiTns of jet emission (e.g. Abdo et al. 2009, 2010ab). FR Is are 
supposed to be the parent population of BL Lac objects, however, 
the Lorentz factors required by a jet model are much lower than 
typical values found in models of BL Lac objects (see, e.g., Abdo 



et al. 2010b). Therefore, the radiation observed in FR Is and BL 
Lacs must have a different origin, which adds some complexity to 
the jet model for FR Is. At least two FR I galaxies reported in Abdo 
et al. (2010b), M87 and Centaurus A, are detected because of their 
proximity rather than a small inclination angle of their jets and, 
then, they are interesting targets for searching for the y-ray emis- 
sion from accretion flows. 

In Niedzwiecki et al. (2012a) we roughly compared prelim- 
inary results of our model with the X/y-ray observational data of 
Centaurus A from INTEGRAL and Fermi-LfiJ'. We found that the 
ADAF model, which matches the X-ray emission in this object, 
predicts the y-ray emission significantly weaker than measured by 
Eermi if the value of the spin parameter a is small. We should 
emend, however, that this conclusion is valid only for models as- 
suming a significant y-ray emission only from the central ~ 10/?g, 
such as our models T and H. Regardless of the values of a and 6, 
our model N with s = 2.6 predicts the absorbed Lo.i-ioGcV approxi- 
mately consistent with 1 .3 x lO'*' erg/s measured in Cen A by Fermi 
(Abdo et al. 2010a). Also regardless of the value of a, our models 
with 5 = 10"' predict the 2-10 keV flux as well as the X-ray spec- 
tral index consistent with that observed in Cen A; for <5 ^ 0. 1 the 
model overpredicts the flux and hardness of the X-ray radiation. 
In the above we assumed the black hole mass of 2 x 10^ Mq (e.g. 
Marconi et al. 2001), which is also within the range allowed by the 
recent measurement of Gnerucci et al. (2011). 

The central core of the accretion system in M87 has been con- 
sidered as the y-ray emitting region e.g. by Neronov & Aharo- 
nian (2007) in their model with particle acceleration in the black 
hole magnetosphere. The accretion rate assessed from the high- 
resolution observations of the nucleus of M87 by Chandra, and 
used to model the multiwavelength spectrum of the nucleus by 
emission from an ADAF, by Di Matteo et al. (2003; note that their 
definition of m differs from ours by a factor of 10), ^ 0. IMg/year, 
corresponds to m =^ 0.01 for a black hole with M = 3 x 10' Mq. 
At such m the central region should be transparent to y-ray pho- 
tons. Then, we compare the unabsorbed luminosities from our non- 
thermal models with s = 2.2 (approximately consistent with the 
slope of the Eermi data above 200 MeV; Abdo et al. 2009) and 
M = 3 X 10' Mg with the y-ray measurements of M87. We find that 
the luminosity derived in the 0.2-10 GeV range by Eermi (Abdo 
et al. 2009), and above 100 GeV by ground-based telescopes (e.g. 
Aleksic et al. 2012) in the low state of M87, can be reproduced by 
our model with a = 0.998 and ^ 0.14Mo/year (for 6^,^ = 40°). The 
required accretion rate, larger by 40 per cent than the face value of 
the estimate in Di Matteo et al. (2003), is allowed by the precision 
of the estimation of the accretion rate using Bondi accretion theory. 
We conclude that hadronic processes in ADAF can contribute sig- 
nificantly to the y-ray emission observed in M87 in the low state, 
or even explain these observations entirely. However, such a model 
requires that all the available material forms the innermost flow, i.e. 
it does not allow a strong reduction the accretion rate in the central 
region by outflows (assumed in models of Sgr A*, see below). 

7.2 Sgr A* 

The value of m in the central region is the major issue for applica- 
tions of ADAF models to the supermassive black hole in the Galac- 
tic Center. The Bondi accretion rate corresponds to m^ ^ 10"' for 
M ^ 3 X 10* Mq and early models used such m to explain the broad- 
band spectra of Sgr A* (e.g. Narayan et al. 1998). The measurement 
of the millimetre/submillimetre polarization of Sgr A* is often as- 
sumed to limit the accretion rate to m < 3 x 10"' (e.g. Marrone et al. 
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2007), but counterarguments are presented, e.g., by Moscibrodzka, 
Das & Czemy (2006) and Ballantyne, Ozel & Psaltis (2007). An 
updated ADAF model to Sgr A* , with fn — 10"^ and strong electron 
heating, is presented in Yuan, Quataert & Narayan (2003). An al- 
ternative to ADAF model, proposed by Moscibrodzka et al. (2006) 
and recently applied by Okuda & Molteni (2012), with the low an- 
gular momentum flow, assumes m ^ 6 x lO""*. The low angular 
momentum flow has a similar density in the innermost region as 
an ADAF with the same m, but much smaller proton temperature 
< 4 X 10" K. Then, the magnitude of hadronic processes could be 
used to distinguish the two classes of models; however, the low T 
may be an artificial effect resulting from the neglect of viscosity in 
the former (low angular momentum) class. 

The CGROIEGKET source, 3EG J1746-2851, was initially 
considered as a possible y-ray counteipart of Sgr A* ; Narayan et al. 
(1998) and OM03 found that their ADAF models underpredicted 
the y-ray luminosity implied by the EGRET measurement. How- 
ever, improved analyses (e.g. Pohl 2005) subsequently indicated 
that 3EG J1746-2851 is displaced from the exact Galactic Center 
and excluded Sgr A*, as well as the TeV source observed with 
HESS which may be directly related with Sgr A*, as its possible 
counterparts. 

Recently, Chemyakova et al. (2011) analysed the Fermi-LJKI 
observations of the Galactic Center and combined them with the 
HESS observational data of the point-like source (Aharonian et al. 
2009). The spectrum of the source seen in the MeV/GeV band by 
Fermi is consistent with a ;r''-decay spectrum. The HESS data in- 
dicate flattening of the spectrum above 100 GeV, suggesting that 
at the highest energies a different spectral component dominates, 
with rather hard spectrum, F ^ 2.2. In the EFg plot, the normal- 
ization of the HESS source is an order of magnitude larger than 
the quiescent X-ray emission of Sgr A*. To explain these observa- 
tions, Chernyakova et al. (2011) discuss a model, following previ- 
ous works (e.g. Atoyan & Dermer 2004), with protons accelerated 
in the accretion flow and then diffusing outwards to interact with 
dense gas at distances of ~ 1 pc. The emission at energies below 
100 GeV may be explained by interactions of protons injected into 
the interstellar medium during a strong flare of Sgr A* that occurred 
300 years ago, as protons generating photons with such energies are 
still diffusively trapped in the y-ray production region. On the other 
hand, most of the higher energy protons have already escaped and 
hence an additional, persistent injection of high-energy protons has 
to be assumed to account for the HESS observations, with the re- 
quired rate of 2 x 10^' erg/s implying a very high efficiency of the 
conversion of the accreting rest mass energy into the proton energy. 

We remark that alternatively the spectral component above 
100 GeV can be explained by proton-proton interactions in an ac- 
cretion flow around a rapidly rotating black hole. Our model with 
a = 0.998, 5 = 2.2 and edge-on viewing direction predict the flux 
consistent with the HESS detection for m = 3 x 10"*. The sce- 
nario with the y-ray emission produced in the accretion flow may 
be tested over the following years, if the fall of the gas cloud into 
the accretion zone of Sgr A* (Gillessen et al. 2012) results in the 
increase of the mass accretion rate. 



7.3 Seyfert galaxies 

Spectral properties of the Seyfert galaxy, NGC 4151, are consistent 
with the model of an inner hot flow surrounded by an outer cold 
disc (e.g. Lubihski et al. 2010), however, its rather large luminos- 
ity, L > O.OlLEdd indicates that a potential y-ray emission from 
innermost region would be strongly absorbed. The constraint of 



io.i-ioGcv/ii4-it)5kcV < 0.0025 derived for this object in Ackermaun 
et al. (2012) can still give interesting information if the spectral 
model constrains the parameters of the inner hot/outer cold accre- 
tion system (in particular, the distance of transition between the two 
modes of accretion). The constraint of Lo.i-ioGcv/i'i4-i95keV < 0.1, 
or even < 0.01 for some objects, found for other Seyfert galaxies 
by Ackermann et al. (2012) is not strongly constraining for the spin 
value, as can be seen in Fig.[8](note that Z-n-igskev is by a factor of 
a few, depending on F^, larger than L2-iokcv)> especially for large 
X-ray luminosities (and hence strong yy absorption) characterising 
most of objects analysed in that paper. 



8 SUMMARY 

We have studied the y-ray emission resulting from proton-proton 
interactions in two-temperature ADAFs. Our model relies on the 
global solutions of the GR hydrodynamical model, same as OM03, 
but we improve their computations by taking into account the rela- 
tivistic transfer effects as well as the yy absorption and by properly 
describing the global Comptonization process. 

We have found that the spin value is reflected in the prop- 
erties of y-ray emission, but the effect is not thrilling. The speed 
of the black hole rotation strongly affects the y-ray emission pro- 
duced within the innermost 1 ORg . If emission from that region dom- 
inates, the observed y-ray radiation depends on the spin parameter 
noticeably; the intrinsic y-ray luminosities of flows around rapidly- 
rotating and non-rotating black holes differ by over a factor of ~ 10. 
However, if the y-ray emitting region extends to larger distances, 
the dependence is reduced. In the most extreme case with protons 
efficiently accelerated to relativistic energies in the whole body of 
the flow, the regions within and beyond 10R„ give comparable con- 
tributions to the total emission, reducing the difference of y-ray 
luminosities between high and low values of a to only a factor of 
~ 2. The radial emissivity profile of y-rays is very uncertain (as the 
acceleration efficiency may change with radius), then, the level of 
the y-ray luminosity cannot be regarded as a very sensitive probe 
of the spin value. Still, it may be possible to assess a slow rotation 
in a low luminosity object by putting an upper limit on the y-ray 
luminosity at the level of < 0.01 of the X-ray luminosity. The pres- 
ence of nonthermal protons may be easily assessed from the y-ray 
spectrum as for a purely thermal plasma the total y-ray emission 
would be observed only at energies lower than 1 GeV. 

We have considered the accretion rate of m = 0.1, for which 
our model gives the bolometric luminosities between ^ 4x lO^^Lndd 
and lO^^Lfidd- Such rii, with the corresponding range of L/Leu, 
seems to be favoured for investigation of the hadronic y-ray emis- 
sion, with the related effects of the space-time metric, because the 
internal y-ray emission is large and its attenuation by yy absorp- 
tion is weak. Flows with L > lO^-^Edd can produce observable 
y-ray radiation only if the emitting region extends out to a rather 
large distance of several tens of Rg. 

We have found the X-ray to y-ray luminosities ratio as a func- 
tion of the black hole spin and the efficiency of the direct heat- 
ing of electrons. The Lo.i-ioGcv/i2-iokcV ratios reaching ~ 0.1, with 
the corresponding levels of the y-ray luminosities which may be 
probed in nearby AGNs at the current sensitivity of Fermi-hAT sur- 
veys, encourage to consider contribution from an accretion flow to 
the y-ray emission observed in low-luminosity objects. We point 
out that such a contribution should be strong at least in Cen A. 

The y-ray luminosity decreases rapidly (oc m^) with decreas- 
ing m. M87 and Sgr A* are the obvious, however possibly also 
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unique, objects for which the y-ray emission from the flow can 
be searched at m « 0.1. Obviously, contribution from other y-ray 
emitting sites should be properly subtracted to establish the lumi- 
nosity of an accretion flow, which may be particularly difficult in 
Sgr A*. The ADAF model with a hard (acceleration index i == 2.2) 
nonthermal proton distribution can explain the y-ray detections of 
both Sgr A* (above 100 GeV) and M87, however, it does not allow 
for strong reduction of the accretion rate by outflows. Nevertheless, 
it seems intriguing that in both nearby, low accretion-rate objects, 
in which the y-ray radiation is not suppressed by yy absorption, 
observations reveal very high energy components, consistent with 
predictions of such a model. In both objects the model requires a 
large black hole spin. 

The luminosity ratio should decrease rather slowly {LyjLx oc 
m"^, as the radiative efficiency of hot flows varies as m"' at small lii; 
cf. Xie & Yuan 2012) with decreasing m. The presence of the y-ray 
signal in low-luminosity AGNs is often considered as the evidence 
for the origin of their radiation in a jet (see, e.g., Takami 201 1). We 
note that such a diagnostic is not valid because accretion flows may 
produce a similarly strong y-ray emission as jets in such objects. 
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